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Abstract 

The decays of ip(2S) into ^ypp, iT pp and npp have been studied with the CLEO-c detector using 
a sample of 24.5 million ip(2S) events obtained from e + e~ annihilations at s/s = 3686 MeV. The 
data show evidence for the excitation of several iV* resonances in pir° and prj channels in ir°pp and 
npp decays, and f2 states in jpp decay. Branching fractions for decays of ip(2S) to ^pp, 7T°pp and 
npp have been determined. No evidence for pp threshold enhancements was found in the reactions 
tp(2S) — > Xpp, where X = 7,7r°,?7. We do, however, find confirming evidence for a pp threshold 
enhancement in J/ip — > "fpp as previously reported by BES. 



* Now at: Pacific Northwest National Laboratory, Richland, WA 99352 
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I. INTRODUCTION 



There is long-standing interest in 6-quark dibaryons and 3 quark - 3 antiquark "baryo- 
nium" states which are permitted in QCD, and may possibly exist. Of particular interest is 
a possible bound state of a proton and an antiproton. The pp state, sometimes called "pro- 
tonium", was searched for in many experiments, but no credible evidence was found [1, 2]. 
Interest was revived in 2002 by two reports by the Belle Collaboration of threshold en- 
hancements in M{pp) in the decays B ± — > K ± pp [3] and B° — > D*°pp [4]. These reports 
were followed by a BES report of threshold enhancement in the decay J/ip — > ^pp [5]. 
Subsequently, there have been reports of threshold enhancements and studies by Belle in 
B + — > n + pp, B° — > K°pp, and B + — > K* + pp [6]; by BaBar in B + — > K + pp [7]; and 
B° ->■ pp + (D°, D*°, D~n + , or D*~n + ) [8]; and most recently by Belle in B + ->• K + pp 
and B + — > n + pp [9]. Many theoretical explanations, cusp effects, final state interactions, 
quark fragmentation, and real bound states of quarks and gluons, have been suggested for 
these threshold enhancements [10]. 

If the enhancement reported by BES in the decay J/ip — > ^pp [5] is due to a threshold 
resonance, it is reasonable to expect that evidence for it may be found also in ip(2S) — > ^pp. 
Further insight into its nature may be provided by the study of the reactions ip(2S) — > n°pp 
and rjpp. 

II. EVENT SELECTION 

In this paper we report on studies of these reactions observed in the CLEO-c detector in a 
data sample of 24.5 million ip(2S) events obtained by e + e~ annihilations at y/s =3.686 GeV 
at the Cornell Electron Storage Ring, CESR. In addition, we use 20.7 pb _1 of off-resonance 
data taken at ^/s = 3.67 GeV. 

The CLEO-c detector, described in detail elsewhere [11], has a solid angle coverage of 93% 
for charged and neutral particles. The charged particle tracking and identification system 
operates in a 1.0 T solenoidal magnetic field, and consists of an inner drift chamber, a 
central drift chamber, and a ring-imaging Cherenkov (RICH) detector. It has a momentum 
resolution of ~0.6% at momenta of ~1 GeV/c. The Csl elecromagnetic calorimeter has a 
photon energy resolution of ~2.2% for = 1 GeV and ~5% at 100 MeV. 

Photons and charged particles with | cos#| < 0.8 were accepted in the detector, where 6 
is the polar angle with respect to the incoming positron beam. For the modes involving the 
direct decays of the ip(2S), exactly two oppositely charged tracks were required in candidate 
events. A photon candidate was defined as a shower which does not match a track within 
100 mrad, is not in one of the few cells of the electromagnetic calorimeter known to be noisy, 
has the transverse distribution of energy consistent with an electromagnetic shower, and has 
an energy more than 30 MeV. For 'ypp the number of showers was required to be > 1, and 
it was required that the most energetic shower (the signal photon candidate) does not make 
a 7T° or i] with any other shower with a pull mass < 3o\ For -n°pp and r/pp the number of 
showers was required to be > 2. 

To identify charged tracks as protons and antiprotons, the energy loss in the drift cham- 
bers (dE/dx) and RICH information was used. For tracks of momentum less than 0.9 
GeV/c, only dE/dx information is used. To utilize dE/dx information, for each particle 
hypothesis, X = it, K, p or p, we calculate = \{dE / dx) mcas — (dE/dx) prc d\/crx, 

where (dE / dx) meas is the measured value of dE/dx, (dE / dx) prc d is the predicted value for 
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hypothesis X , and ax is the standard deviation of the measurements for hypothesis X . We 
cut on both the deviation of the measured dE/dx from a given particle hypothesis, \ 
and the difference in ■^ dE / dx between two particle hypotheses, Ax^y** = Xx^ dX ~ Xy^*- 
For higher momentum tracks, we use a combined log-likelihood variable. For example, to 
differentiate between proton and pion we construct 

A£ Pi7r = ( X f 1*°)* - (xf^y + 2 x {Lf CH - L« ICH ) 

where are the log-likelihoods obtained from the RICH subdetector. We use RICH 

information if the track has |cos#| < 0.8 and the track has valid RICH information for 
at least one hypothesis (pion or proton), and at least three photons consistent with that 
hypothesis were recorded in the RICH. 

We consider three different momentum regions for charged tracks. 

• p < 0.9 GeV/c: In this momentum region only dE/dx information for the tracks is 
available, and it is required that it be within 3a p of the proton hypothesis, and must 

be more "proton-like" than "pion-like" or "kaon-like", i.e. \xp E ^ dx \ < 3, Axn,p > 0, 
and Axk, p > 0. 

• 0.9 GeV/c < p < 1.15 GeV/c: In this momentum region, although we are above 
the threshold for a proton to emit Cherenkov radiation in the RICH, the probability 
that it will do so is still low. Therefore, if RICH information is available, we require 
that the track be more "proton-like" in the combined log-likelihood variable, i.e., 
A£ p 7r < 0. If RICH information is not available, we again require that \xp ^ 1^3, 
and additionally require a 5a difference between the proton hypothesis and the pion 
and kaon hypotheses, i.e., Ax n , P > 5, and Axk, p > 5, in order to reduce the number 
of other particles which pass these cuts. 

• p > 1.15 GeV/c: In this momentum region, dE/dx alone no longer provides useful 
information for proton identification. We require that RICH information be available, 
and that A£ Pj7r < 0. 

We require one of the charged tracks to be identified as a proton or antiproton and 
assume the other track to be its antiparticle as required by baryon conservation, and we 
require the proton and antiproton to come from a common vertex, with kinematic fit yielding 

v 2 - < 20 

App vertex ^ 

Finally, in order to select the events for the channels of interest: 

• For selection of ^>(2S) — > ^pp events we require xit/degrees of freedom (d.o.f.) < 5 
for the four-momentum conservation constrained fit to p, p and the signal photon 
candidate. 

• For selection of tp(2S) — > ir°(ri)pp events, we first require that only one 7r°(?7) be made 
by any two photons and the pull mass be within 3a. Then we require xlt/d-o.f. < 5 
for the four- momentum conservation constrained fit to p, p, and 7r°(r/). We remove 
the events corresponding to ip{2S) — > 7r°(r))J/ij) by rejecting candidates for which 
M{pp) = M(J/ip) ±20 MeV/c 2 . Figure 1 shows the distribution of M^j) before and 
after the selection of 7r° and r\ described above. 
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• For selection of ip(2S) —> n + n J/if), J/4> —> •ypp events the additional event selection 
requirements are described in Sec. VIII. 

The values of the \ 2 cu ts for the fits were selected based on the comparison of the data 
and the phase space distributions for the individual decays obtained from Monte Carlo (MC) 
simulations. 

3850610-001 




M (yy) (GeV/c 2 ) 

FIG. 1. Distribution of M{^yy) in ip(2S) — > ir°(r])pp. The unshaded histogram shows the M("fj) 
distribution before the selection of ir° and 77, and the shaded histogram shows it after the n° and 
rj selection described in the text. 



III. MONTE CARLO STUDIES 

In order to verify the event selection criteria and determine efficiencies, 50,000 phase space 
MC events were generated for each decay channel analyzed. As an example, for ip(2S) — > ypp 
events with M(pp) <2.85 GeV/c 2 the contribution of each step of event selection is presented 
in Table I. The overall phase space efficiency is (27.7 ± 0.2)%. The corresponding efficiency 
for 7p(2S) n°pp is (26.9 ± 0.2)%, and for ip(2S) rjpp it is (27.7 ± 0.2)%. 

We also use ip{2S) "generic" MC events with the available statistics of about five times 
data events (~118 million events). The generic if>(2S) MC sample is generated using the 
available branching fractions for the i/j(2S), Xcj, J/i>, and r\ c decays, with unmeasured 
decay modes simulated by JETSET [12]. We have tested the event selection using a generic 
MC sample. We apply the same event selection to these MC events, extract the different 
branching fractions, and compare them to the branching fractions which were input in 
creating the generic MC sample. The agreement between the input and output branching 
fractions for ip(2S) —> •ypp, 7i°pp, and r/pp, is found to be within (2.4 ± 3.9)%, (1.0 ± 1.0)%, 
and (3.0 ± 4.0)%, respectively. 

IV. OVERVIEW OF ip(2S) DECAYS 

Figure 2 shows Dalitz plots for the data for the three decays. All three Dalitz plots 
show event populations which are far from uniform, as would be expected for pure phase 
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TABLE I. Efficiencies of the individual event selection criteria for the decay ip(2S) — > jpp based 
on phase space MC simulation. 



Selection requirement 


tmciency (to J 


Charged track and photon selection 




N c h = 2, net charge = 0, iV 7 >1 


74.5 


Signal photon does not make 




7T° with any other shower 


97.0 


Vertex fit, constrained fit 




Xvertex < 2U > Xfit < 




Proton-antiproton identification 


98.2 


M(pp) <2.85 GeV/c 2 


64.3 


All p, p and most energetic photon 




are in the barrel ( cos# < 0.8) 


72.0 



Total 27.7 
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FIG. 2. (a) Dalitz plots for the data: (a) M 2 (pj) versus M 2 (pj) for Y>(2S) ->• 1P p- (b) M 2 (pvr°) 
versus M 2 (p7r°) for ijj(2S) -)• 7r°pp; (c) M 2 {pij) versus M 2 (pr]) for ^(25*) ->■ 



space decays, and suggest contribution by intermediate excited nucleon states, N*, and 
mesons. Since the branching fractions for A^* decays to Nir and Nrj are generally much 
larger than those for decays to Nj, we expect excitation of N* states in i/)(2S) — > n pp and 
ip(2S) — > rjpp. Similarly, we expect excitation of intermediate meson states which decay into 
pp, fj states in ip(2S) — > ^pp and rjpp, and aj states in ip{2S) — > n°pp. 

A caveat about the intermediate states is in order. The intermediate N* and / 2 states 
which we use in our analysis tend to have masses > 1.5 GeV/c 2 . Unfortunately, the 
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existence of most high mass N* resonances is either uncertain or poorly established, and 
their masses and widths have large uncertainties, so much so that the 2008 Particle Data 
review (PDG08) [13] omits many of them from its summary table. Similar uncertainties exist 
for meson states with masses > 1.5 GeV/c 2 . Therefore our identification of an observed 
resonance with a known resonance is necessarily tentative. 

Our data lack the statistics to make a full partial wave analysis. Instead, we analyse the 
projections of the Dalitz plots of invariant mass distributions for M(p(7, 7r°, T])) and M(pp). 
Throughout this paper, charge conjugate states and their contributions are implied. 

We fit the invariant mass distributions with contributions from phase space and the min- 
imum number of resonances required to obtain good fits. The resonances are parametrized 
in terms of relativistic Breit-Wigner functions with mass dependent widths and include the 
Blatt-Weisskopf penetration factors [13, see p. 772]. We note that peak positions and widths 
in the relativistic fits can be substantially different from those for the simple Breit-Wigner 
function, particularly for large widths and proximity to thresholds [13]. 

In order to take proper account of intermediate states and possible reflections in the Dalitz 
plots we analyze the data in the full range of pp invariant mass, from threshold to 3.6 GeV/c 2 
. In this important respect the present analysis differs from the BES analyses [5, 16, 17]. 

In the following Sections, V, VI and VII we discuss the decays ip(2S) —> jpp, Tr°pp and 
f]pp, respectively. In Sec. VIII we present the results of the analysis of our limited statistics 
sample of J /if) — > •ypp events. 



V. THE DECAY i/>(2S) ->• -ypp 

For the search for a threshold enhancement in ip(2S) — > jpp, and for the measurement 
of the inclusive branching fraction B(tp(2S) — > ^pp) we limit ourselves to M(pp) < 2.85 
GeV/c 2 , below the r\ c mass. As a check on our analysis technique, we use our data for 
M(pp) > 3.15 GeV/c 2 to calculate the B(xc,j PP) branching fractions and compare them 
to recent measurements. 



3850610-003 




(a) 



UL ! ! ! LDJ L- 




20 40 60 80 100 20 40 60 80 100 



' pp vertex 



FIG. 3. Distributions of the x 2 °f four-momentum conservation constrained fits for ip(2S) — > jpp: 
(a) % 2 of vertex fit, (b) x 2 °f the four-momentum conservation fit. Points correspond to the data 
and the shaded histograms correspond to the phase space MC simulation. Dashed lines indicate 
the cut values used. 
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FIG. 4. Distributions of E/pc for ip(2S) — > ^ypp for M(pp) < 2.85 GeV/c 2 : (a) protons, (b) 
antiprotons. Because of annihilations E/p for antiprotons extends over a much larger range than 
for protons. Points correspond to the data and the shaded histograms to the phase space MC 
simulation. 




FIG. 5. Distributions of the momenta of charged particles for ip(2S) — > "fpp: (a) protons, (b) 
antiprotons. Points correspond to the data and the shaded histograms to the phase space MC 
simulation. 



In Fig. 3 we show the \ 2 distributions for the data and the phase space MC events 
for vertex fit and four-momentum conservation constrained fit to the proton, antiproton 
and most energetic shower in the event. All other selection criteria have been applied. 
Comparison of these distributions suggests the cut values Xpp vertex < 20 and xlt < 20. 

In Fig. 4 we compare E/pc distributions for protons and antiprotons in data and the 
phase space ip(2S) — > 'ypp MC simulation, where E is the energy determined from the 
calorimeter and p is the momentum determined from track reconstruction. The distributions 
of protons and antiprotons are different, because antiprotons annihilate in the material of 
the electromagnetic calorimeter. However, for both protons and antiprotons, the data and 
the phase space MC distributions show good agreement. 

In Fig. 5 proton and antiproton momentum distributions for the data and the phase space 
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TABLE II. Results for B{x c j -> Pp) for Xcj states. The PDG08 values for B(ip(2S) -»■ 7^ c j) have 
been used to obtain the results for B(x c j — > PP)- 







XcO 


Xcl 


Xc2 


Mass (MeV/c 2 ) 




3412.8±1.0 


3512.5±0.4 


3555.0±1.0 


N (events) 




236.0 ± 18.4 


79.0 ± 10.7 


62.5 ± 9.8 


Efficiency in % 




39.5 


41.4 


37.9 


B(if>(2S) -> 7 Xcj) x B( Xc j -> pp) x 


10 6 22.0 ± 1.7 


7.05 ± 0.96 


6.07 ± 0.95 


B(XcJ -»■ PP) x 10 5 


(this analysis) 


23.4 ± 2.1 


8.0 ± 1.1 


7.3 ± 1.2 


B(XcJ ->■ PP) x 10 5 


(CLEO [14]) 


25.7 ± 1.5 ± 2.0 9.0 ± 0.8 ± 0.6 7.7 ± 0.8 ± 0.6 


B(Xcj -> pp) x 10 5 


(PDG08 [13]) 


21.4 ± 1.9 


6.6 ± 0.5 


6.7 ± 0.5 



ip(2S) — > 7pp MC simulation are shown. The agreement between MC simulation and data 
momentum distributions is not good, and may indicate the effect of intermediate resonances. 

In Fig. 6 we present the M(pp) invariant mass distributions for the data and the ip{2S) 
generic MC simulation, which includes the excitation of Xcj and r) c , but not the ISR generated 
J/i/j. All event selection criteria have been applied. The generic MC events are normalized 
to the number of ip{2S) events in the data for a qualitative comparison. 
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FIG. 6. Distribution of M{pp) for iJj(2S) "iPP- All event selection criteria have been applied. The 
points represent data and the shaded histogram is the generic MC distribution, which is normalized 
to the 24.5 million tp(2S) events in data. The dashed line indicates the M(pp) = 2.85 GeV/c 2 limit 
of the range used for ip(2S) — > jpp branching fraction calculations. 

In Fig. 6 we note that there is an excess of events in the data over the generic MC 
simulation in the region M(pp) < 2.3 GeV/c 2 . The generic MC simulation has no input for 
possible resonances in this mass region. 
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FIG. 7. Fit of the photon energy distribution for ip(2S) — > jpp. Events in the figure are the same 
as in Fig. 6. The fit has x 2 /o!-o-.f. = 100/98. 



The excitation of the Xcj states, shown in Fig. 6 gives us an opportunity to further test 
the appropriateness of our event selection. 

To extract Xcj branching ratios we fit the photon energy E 1 distribution shown in Fig. 7 
with Breit-Wigner functions convolved with Crystal Ball line shape [15], and a second-order 
polynomial background. The fit results for photon energies agree with those expected for the 
Xcj resonances within ±2 MeV. As seen in Table II, our calculated values of B(xcJ ~^ Pp) 
agree within errors with both the results of a recent CLEO measurement [14], which use the 
same data, and PDG08 [13]. These measurements are intended as checks on the analysis 
technique only and not as new measurements, and no systematic errors are included. 

In order to explore the intermediate state resonances which are excited in the reaction 
ip(2S) — > '-/pp we study several presentations of the data. These include the py versus jry 
Dalitz plot, the M(pp) and M(pj) projections, and the distributions of cosO, where is 
the angle between the proton and antiproton in the rest frame of the photon-proton system. 

The three panels in Fig. 8 show M 2 (jry) versus M 2 (jry) Dalitz plots respectively for 
phase space MC simulation, data, and MC simulation with the intermediate resonances as 
described below. Figure 9 shows phase space MC distributions superimposed on the data 
for M(pp), M(p7,p7) and cos©. (In the last two plots events have been included for both 
P7 and pj, resulting in double counting). 

In the distributions shown in Figs. 8 and 9, it is clear that pure phase space distributions 
fail to describe the data. Significant contributions by intermediate states are required. We 
have made MC studies of the contributions that various known scalar and tensor meson reso- 
nances would make to these distributions. We find that the best candidates are /2(1950) and 
/2(2150) with parameters given in PDG08. We determine MC shapes of the contributions 
that / 2 (1950), / 2 (2150), and phase space make to the distributions for M(pp), M(jry) and 
cos0, and determine their relative magnitudes by fitting the data distributions. Using the 
PDG08 [13] values for masses and widths for the two resonances, the best fit is obtained with 
relative fractions listed in Table III. The corresponding MC determined efficiencies, which 
are found to be insensitive to the uncertainties in masses and widths of the resonances, are 
also listed in the table. The overall efficiency of the admixture of the resonances and phase 
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FIG. 8. Dalitz plots of M 2 (p 7 ) versus M 2 (p r y) for ijj(2S) — > 'ypp: (a) the phase space MC simulation; 
(b) the data; (c) the sum of three MC plots for /2(1950), /2(2150) and phase space. 



385061 0-009 




M (pp) (GeV/c 2 ) M (py) (GeV/c 2 ) cos 



FIG. 9. Distributions in data (points) compared to the phase space MC distributions (solid lines): 
(a) the M{pp) invariant mass distributions; (b) the M{p^f) distributions; (c) the cos distributions. 
Phase space normalization is arbitrary in all plots. 




FIG. 10. Distributions in data (points) compared with the sum of the MC distributions in the 
proportions given in Table III (solid lines): (a) the M{pp) invariant mass distributions; (b) the 
M(jprj) distributions; (c) the cos distributions. The individual contributions are: /2(1950) shown 
with the dotted line marked 1, /2(2150) shown with the dotted-dashed line marked 2, and phase 
space shown with the dashed line marked 3. 
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TABLE III. Fractions and efficiencies for the intermediate resonances and the phase space contri- 
bution for the best fits for the reaction ip(2S) — > jpp. 



State 


M (MeV/c 2 ) T (MeV/c 2 ) Fraction (%) 


e 


/ 2 (1950) 


1944±12 472±18 32±5 


0.375 


/ 2 (2150) 


2156±11 167±30 21±5 


0.410 


Phase space 


47±6 


0.277 



space is 

(e) = 0.336 ±0.008. (1) 

As shown in Fig. 10 good fits to all three distributions are obtained with the above 
admixtures, their respective x 2 /d.o.f. being 15/24 for M(pp), 33/35 for M(jry), and 23/20 
for cos 0. The resulting Dalitz plot, shown in Fig. 8(c), is also in good qualitative agreement 
with that for the data. No evidence is found in the data for a narrow resonance R with 
T R < 40 MeV/c 2 anywhere in the region M(pp) = 2200-2800 MeV/c 2 . The 90% confidence 
level upper limit is B(if>(2S) ->■ jR) x B{R ->■ pp) < 2 x 10~ 6 . 



A. Determination of B(tp(2S) — > -ypp) 

In the region M(pp) < 2.85 GeV/c 2 , we obtain A" = 407 ± 20 ip{2S) — > ^pp candidate 
events. We evaluate the background due to ip{2S) — > n°pp, in which one photon from 
the 7T° decay is lost, as A" bkg (7r ) = 38 ± 2. In addition, by analyzing the continuum data 
at y/s = 3.76 GeV we determine that the luminosity- normalized continuum background 
contribution is N cont = 26 ±8 counts. With the relative contributions of / 2 (1950), / 2 (2150), 
and phase space as in Table III, and the effective overall efficiency (Eq. (1)) we get 

Bms) ^ m) = N ~ Ncoat 

{€) X 1\^2S) 

= (4.18 ± 0.26(stat)) x 10" 5 . (2) 

The individual product branching fractions are 

B(i/>(2S) 7/ 2 (1950)) x B(/ 2 (1950) ^ pp) 

= (1.2±0.2(stat)) x 10~ 5 , (3) 



B{^{2S) -)■ 7/2(2150)) x 8(/ 2 (2150) -)• pp) 

= (0.72 ± 0.18(stat)) x lO" 5 . (4) 

Estimates of systematic errors are provided in Sec. IX. 

Our result for B(tp(2S) — > •ypp) differs by 2a from the PDG08 result based on the BES 
measurement [16], B(ip(2S) — > ^pp) = (2.9 ±0.6) x 10 -5 in which no account of intermediate 
resonances was taken. The results in Eqs. (3) and (4) represent the first measurements of 
these product branching fractions. 
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FIG. 11. (a) MC-determined efficiency as a function of AM = M{pp) — 2m p for ip(2S) — > jR t h r , 
Rthr — > PP- (b) The solid curve is the shape of the AM distribution from MC simulation for the 
admixtures shown in Table III, and the dashed curve is the shape of the AM distribution from 
MC simulation for phase space alone. Relative normalizations are arbitrary. 



3850610-012 




0.05 0.10 0.15 0.20 0.25 0.30 



AM = M(pp) - 2m p (GeV/c 2 ) 



FIG. 12. Fits of the AM = M(pp) - 2m p distribution for ip(2S) -fpp decay. The dotted line 
is the sum of the resonance and phase space contributions according to Table III, the dashed line 
shows the fitted threshold resonance contribution. The solid line is the sum of all contributions. 



B. Search for Threshold Enhancement in rp(2S) — > ~/pp 

Figure 10(a) shows that a good fit to the M{pp) spectrum is obtained with the sum 
of contributions from / 2 (1950), / 2 (2150), and phase space, with xVd.o.f. = 15/24. No 
threshold resonance seems to be needed. However, to reach a quantitative conclusion we 
study in detail the pp threshold region, AM = M(pp) — 2m p = — 300 MeV/c 2 . To do 
so, we evaluate the contributions of the / 2 resonances and phase space in this region, and 
determine the efficiency for ip{2S) — > 7-fiW, -Rthr — > PP m this region. The results are shown 
in Fig. 11. 
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The AM = M(pp) — 2m p event distribution is shown in Fig. 12. A visual inspection of 
the distributions shows that there is no evidence for a statistically significant enhancement 
at the threshold, AM = 0. In fact, a straight line fit to the data gives x 2 /d.o.i. = 52/58. 
However, we must consider the contributions due to the / 2 resonances and phase space as 
has been determined in Table III, and as shown in Fig. 11(b), and the efficiency (e) = 55.8% 
in the threshold region. Figure 12 shows the best fit obtained using these contributions 
plus a Breit-Wigner threshold resonance with the parameters obtained by BES [5], namely 
M{pp) = 1859 MeV/c 2 , and T = 20 MeV/c 2 . The fit has x 2 /d.o.f. = 53/58, and includes 
the best fit threshold resonance R t hr with 9±g° counts. This leads to 

BmS) 7^thr) x B(R thr pp) = 
= (0.66+^) x 1CT 6 or < 1.6 x 1(T 6 , 90% CL, (5) 

where CL means confidence level. This is more than a factor three more restrictive than the 
current best limit [16]. 



VI. THE DECAY ^(25) -»■ 7r°pp 

Our analysis of ip(2S) — > ir°pp follows the same steps as described in Sec. IV for tp(2S) — >■ 
'-/pp. Figure 13 shows the three Dalitz plots respectively for (a) phase space MC simulation, 
(b) data, and (c) MC simulation with the resonances described below. The phase space and 
data Dalitz plots differ dramatically, and the MC plot with the resonances described below 
is in impressive agreement with the data. Figure. 14 shows the projected distributions for 
M(j07r°), M(pp), and cos 0, the polar angle of p in the rest frame of n°p. It is clear that the 
pure phase space distributions do not reproduce the data in either the Dalitz plots or the 
projected distributions. Figure 15 shows the same three distributions with good quality fits 
based on resonance shapes determined from MC simulations, as described below. 

The N* intermediate states in tp(2S) — > ir°pp are most clearly visible in the M(pn°) 
distribution of Fig. 14(a), with enhancements near M(pn°) m 1400 MeV/c 2 and M(p7r°) m 
2300 MeV/c 2 . Similarly, the meson intermediate states are most clearly visible in the M(pp) 
distribution of Fig. 14(b), with enhancements near M(pp) ^2100 MeV/c 2 amd M(pp) w 
2900 MeV/c 2 . The enhancement at MQott ) w 1400 MeV/c 2 can be identified with the 
well known A*(1440), which we call N£, and the enhancement at M(pp) fa 2100 MeV/c 2 , 
which we call Ri, can be identified with the known resonance /o(2100) [13]. The large 
enhancements in MQtwr ) at 2300 MeV/c 2 , which we call N*, and in M(pp) at 2900 MeV/c 2 , 
which we call R 2 , can not be identified with known N* and f 2 resonances, and we have to 
take an empirical approach for them. 

Because the mass and width of /o(2100) are well defined, in all subsequent analysis we 
keep them fixed to their PDG08 values. To determine the optimum values for the masses 
and widths of the N£, iV| and R 2 resonances following procedure was used. 

Because the Dalitz plot projections contain reflections, the projections can not be fitted 
with simple Breit-Wigner resonances. Instead, MC distributions have to be generated for 
individual resonances with assumed masses and widths, and their optimum values have to 
be determined by fitting the data distributions with the MC generated distributions. Our 
procedure takes account of reflections, but does not include taking account of any possible 
interferences between resonances. 
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FIG. 13. Dalitz plots of M 2 (p7r°) versus M 2 (p7r°) for the channel i[i(2S) — > 7r°pp: (a) the phase 
space MC simulation; (b) the data; (c) the sum of four MC plots for i?i(2100), i? 2 (2900), iY\*(1440) 
and iV| (2300). 




FIG. 14. Distributions in data (points) compared to the phase space MC distributions: (a) the 
M(p-K°) distributions, (b) the M(pp) invariant mass distributions, (c) the cos distributions. 
Normalization is arbitrary in all plots. 




FIG. 15. Distributions in data (points) compared to the sum of the MC distributions in the pro- 
portions given in Table IV (solid lines): (a) the M(p7r°) distributions; (b) the M{pp) distributions; 
(c) the cos© distributions. The individual contributions are: i?i(2100) shown with the dashed line 
marked 1, N* (1440) shown with the dotted line marked 2, iV|(2300) shown with the dotted-dashed 
line marked 3, i?2(2900) shown with the dotted line marked 4. 
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TABLE IV. Fractions and efficiencies for the intermediate resonances for the reaction ip(2S) — > 
■K°pp. 

Resonance M (MeV/c 2 ) T (MeV/c 2 ) Fraction (%) e 



^(1440) 


1400±25 


220±20 


50±4 


0.241 


A 2 *(2300) 


2300±25 


300±30 


28±4 


0.276 


#1(2100) 


2103±8 


209±19 


8±3 


0.275 


#2 (2900) 


2900±20 


250±25 


24±4 


0.241 



We first fit the M(pn°) distribution with only and N% resonances and determine the 
best fit values for their masses and widths by iterating each in 5 MeV/c 2 steps. We then fit 
the M(pp) distribution with just the above and A^ resonances. We find that the M(pp) 
distribution is fitted poorly, with x 2 /d.o.f.=65/37, and the enhancements at M(pp) « 2100 
MeV/c 2 and M(pp) 2900 MeV/c 2 are not reproduced. We then explicitly introduce fixed 
parameter i?i(2100), and R 2 on whose parameters we iterate to find their best values. As 
expected, the fit to the M(pp) distribution is improved, with x 2 /d.o.f. =44/33. We go back 
to the M(p7r°) distribution to determine the effect of including Ri and R 2 . It is found that 
their contribution is structureless in the M(pir°) distribution, and it does not affect the best 
fit parameters of N£ and N%. 

In Table IV final resonance parameters of N±, N£, R± and R 2 are listed. The errors in the 
masses and widths are those which change the likelihood of fits by two units. The efficiencies 
are as determined by MC simulations. The relative fractions are determined by the final fit 
to the M{pp) distribution. 

The fits obtained for M(p7r°), M(pp) and cos 6 distributions with the final set of pa- 
rameters for all four resonances are shown in Fig. 15(a,b,c). The individual resonance 
contributions are shown with dotted and dashed lines. The corresponding composite Dalitz 
plot is shown in Fig. 13(c). It agrees very well with that for data in Fig. 13(b). No evidence 
for a pp threshold enhancement is observed in the M(pp) distribution of Fig. 15(b). 

BES reported the result for ip(2S) — > ir°pp using their sample of 14 million ip(2S) 
events [17]. Since the number of events was almost a factor of two smaller than in the 
present investigation, they were not able to reach any conclusions about intermediate states 
other than to note that there was "indication of some enhancement around 2 GeV/c 2 ." 

A. Determination of B(i^(2S) — > 7r°pp) 

We consider all events in the M(pp) spectrum for M(pp) < 3.6 GeV/c 2 (Fig. 15(b)) 
for determination of the branching fraction B(i(j(2S) — > 7r°pp), except those with M(pp) = 
3.097 ± 0.020 GeV/c 2 , which could come from J/ip production. 

We obtain = 1063 ± 33 ip(2S) — > n°pp candidate events. Using our branching fraction 
of B(tp(2S) ->■ ~fpp = (4.18 ± 0.26) x 10" 5 , and the efficiency determined in Eq. (1) we 
estimate Abkg = 15 ± 1 background counts due to misidentified •ypp events. 

We estimate the non-resonant contribution of 7c°pp production by using the data taken 
at the o&-iIj(2S) resonance energy of = 3.67 GeV. It leads to a luminosity-normalized 
non-resonant contribution in our data, N cont = 105 ± 16. 

The efficiencies determined from the AT*(1440), iV 2 *(2300), i?i(2100) and i? 2 (2900) MC 
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simulations are 24.1%, 27.6%, 27.5%, and 24.1%, respectively. The overall efficiency of the 
admixture of the resonances is (e) = (25.4 ± 0.2)%. 
This yields a branching fraction of 

K>///r>rr\ - 0\ N — N^g — N cont 

nms) -» = (e)x%2S)X 8 B(T o^ 77) 

= (1.54 ± 0.06(stat)) x 10~ 4 . (6) 

This result is in agreement with the PDG08 [13] value of (1.33 ± 0.17) x 10~ 4 , and has a 
factor three smaller error. 

We can also determine the product branching fractions for the ^(1440), 7X^(2300), 
i?i(2100) and i? 2 (2900) resonances by taking account of their respective fractions and effi- 
ciencies given in Table IV. The resulting product branching fractions are 

B(ip(2S) ^pN*(lU0)) x B(N*(1U0) ^ pir°) 

= (8.1 ±0.7(stat)) x 10- 5 , (7) 

B(tp(2S) ^pN*(2320)) x B(N* (2320) ^ pn ) 

= (4.0±0.6(stat)) x 10~ 5 , (8) 

B(ip(2S) -> tt° j R 1 (2100)) x i3(i?i(2100) -> pp) 

= (1.1 ±0.4(stat)) x KT 5 , (9) 

B(ip(2S) ->■ 7r° J R 2 (2900)) x £(#2(2900) ->■ pp) 

= (2.3±0.7(stat)) x 10~ 5 . (10) 

These are the first determinations of these product branching fractions. Estimates of 
systematic errors are provided in Sec. IX. 



VII. THE DECAY ^(2S) -> r] P p 

As shown in Fig. 1, the yield of r\pp (184 counts) is nearly a factor of 6 smaller than for 7i°pp 
(1063 counts). Dalitz plots of M 2 (prj) versus M 2 (prf) are shown in Fig. 16, (a) for pure phase 
space, (b) for the data, and (c) for resonances described below. As in the ip(2S) — > ix Q pp 
case, it is seen that the Dalitz plot for the data (Fig. 16(b)) is completely different from the 
uniformly populated Dalitz plot for phase space MC simulation (Fig. 16(a)). The data are 
clearly dominated by the contribution of intermediate states. As shown in Fig. 16(c) the 
sum of MC simulated contributions of resonances described below reproduces the data very 
well. 

The projected distributions for (a) M(pi]), (b) M(pp), and (c) cos 6, where 6 is the polar 
angle of the p in the rest frame of r/p, are shown in Fig. 17, together with MC generated 
distributions for phase space. As expected, the phase space distributions do not reproduce 
the data distributions. 

In the M(pi]) invariant mass distribution (Fig. 17(a)) no evidence is found for 7X^(1440) 
and N% (2300) resonances seen in the M(pir°) plot of Fig. 15(a), but a large peak is observed 
at M(prj) ~ 1540 MeV/c 2 . This suggests excitation of iV*(1535) J p = \ nucleon resonance 
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FIG. 16. Dalitz plots of M 2 (pr]) versus M 2 (prj) for the channel ip(2S) — > r\pp: (a) the phase space 
MC simulation; (b) the data; (c) the sum of two MC plots for #i(2100) (20%), and AT* (1535) 
(80%). 
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FIG. 17. Distributions in data (points) compared to the phase space MC distributions: (a) the 
M(pr)) distributions; (b) the M(pp) invariant mass distributions; (c) the cos distributions. Nor- 
malization of the curves is arbitrary in all plots. 
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FIG. 18. Distributions in data (points) compared to the sum of the MC distributions in the 
proportions given in Table V (solid lines): (a) the Ai(prj) distributions; (b) the M{pp) distributions; 
(c) the cos0 distributions. The individual contributions are: i?i(2100) shown with the dotted line 
marked 1, iV*(1535) shown with the dashed line marked 2. 
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of PDG08 [13] with M = 1525 - 1545 (MeV/c 2 ) and T = 125 - 175 (MeV/c 2 ), which is 
known to decay into Nrj with a branching fraction of (45-60)%. 

In the M(pp) distribution (Fig. 17(b)) there is a broad enhancement in the 2.7 — 3.0 
GeV/c 2 region which arises mainly as reflection of the iV*(1535) resonance in M(pq). In 
addition there is a narrow enhancement near M(pp) ps 2100 MeV/c 2 reminiscent of the one 
observed in M(pp) from ip(2S) — > n pp decay. 

The optimized masses and widths we obtain for these resonances, and their fractions and 
estimated efficiencies are shown in Table V. It is found that the MC-determined efficiencies 
are insensitive to the uncertainties in masses and widths of the resonances. 

TABLE V. Fractions and efficiencies for the intermediate resonances for the best fits for the reaction 
ip(2S) —7- rjpp. 



Resonance M (MeV/c 2 ) 


T (MeV/c 2 ) 


Fraction (%) 


e 


iV*(1535) 1535±10 


150±25 


80±6 


0.294 


#1(2100) 2103±8 


209±19 


20±6 


0.259 



We fit our data distribution for M(prf) and M(pp) with an admixture of MC simulated 
shapes for these two resonances (the fit result shows that the contribution from the phase 
space MC is consistent with zero). The best fit admixture is: 

(0.20 ± 0.06) x i2x(2100) + (0.80 ± 0.06) x JV*(1535). (11) 

The final results for the M(pi]), M(pp) and cosO distributions are presented in Fig. 18 
with solid lines. The dashed and dotted lines represent the contributions from the individual 
resonances. Good agreement between the data and the fitted distributions is obtained for 
all three distributions, with x 2 /d.o.f. are 35/30 (M(pr))), 30/30 (M(pp)) and 32/20 (cos 6). 

The M 2 (prj) versus M 2 (prj) Dalitz plot in Fig. 16(c) constructed with the MC distributions 
for the above resonance admixture is seen to reproduce very well the Dalitz plot of data in 
Fig. 16(b). 

A. Determination of B(il)(2S) — >• r/pp) 

We consider the entire M(pp) spectrum with M(pp) < 3.077 GeV/c 2 in Fig. 18(b) for 
the determination of B(ip(2S) — > r/pp). 

We obtain N = 184 ± 14 ip(2S) — > rjpp candidate events. We do not find any background 
contribution from feed-down from other decay channels. We estimate the continuum con- 
tribution of r/pp production by using the data taken at y/s = 3.67 GeV. It leads to the 
luminosity-normalized contribution N cont = 30 ± 8 counts. 

The reconstruction efficiencies determined from the i?i(2100) and A^*(1535) MC simula- 
tions are 25.9% and 29.4%, respectively and the overall effective efficiency of the resonances 
admixture is (e) = (28.7 ±0.2)%. 

This yields a branching fraction of 

= (5.6±0.6(stat)) x 10~ 5 . (12) 
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This is in agreement with the PDG08 [13] value of (6.0 ± 1.2) x 10 -5 , and has a factor two 
smaller error. We can determine the product branching ratios for the production of iV*(1535) 
and i?i(2100) resonances by taking account of their respective fractions and efficiencies. We 
obtain 

B(ip(2S) ^pA*(1535)) x £(A*(1535) ->■ pr]) 

= (4.4±0.6(stat)) x 10~ 5 , (13) 

B(tp{2S) ->■ 77i?i(2100) x 6(^(2100) ->■ pp) 

= (1.2±0.4(stat)) x 10" 5 . (14) 

These are the first determinations of these product branching fractions. Estimates of 
systematic errors are provided in Sec. IX. 

VIII. SEARCH FOR pp THRESHOLD ENHANCEMENT IN J /if; -> 1P p 

Although the number of 7r + 7r~ tagged J/ip events in our sample of ip{2S) — > tt + 7t~J/i/j 
is 8.7 million, as compared to the 58 million event J/ip sample of BES II, it is instructive 
to analyze it for the sub-threshold resonance with M(pp) = 1859+i +25 MeV/c 2 reported by 
BES [5]. 

For selection of ip(2S) — > n + 7i~ J/ip, J/ip — > •ypp events, we require first Xpp vertex < 20 
and xit/d-o.f. < 3 of the J/ip mass constrained fit to pp and most energetic shower in the 
event, then Xvertex < 40 and Xg t /d.o.f. < 10 of the four-momentum conservation constrained 
fit to 7r + 7r~ and J/ip. 

In Fig. 19 we show the distribution of pp invariant mass as a function of AM = M(pp) — 
2m p in the extended mass region AM = - 970 MeV/c 2 . We believe that it is essential to 
analyze the data in the extended mass region, because as we have seen for ip{2S) decays, 
higher mass resonances make contributions all the way down to the pp threshold. Further, 
as shown in Fig. 19, a much better estimate of the phase space contribution can be made 
when data in the extended mass region is taken into account. Fig. 19 shows an enhancement 
near pp threshold and a large broad enhancement around AM 200 MeV/c 2 . We therefore 
analyze our data in the extended mass region, AM = — 970 MeV/c 2 , and take account 
of possible resonances other than the one near the pp threshold. Our analysis differs in 
this essential respect from that of BES in which data in the limited region, AM = — 300 
MeV/c 2 , was analyzed, and no account was taken of the enhancement around AM ps 200 
MeV/c 2 . 

We have made an attempt to fit the present AM distribution with a threshold resonance 
plus the complement of resonances and phase space observed in the case of ip(2S) — > 7pp, 
i.e. / 2 (1950) (corresponding to AM ^ 74 MeV/c 2 ), / 2 (2150) (corresponding to AM ^ 224 
MeV/c 2 ), and phase space. No evidence for a contribution due to the / 2 (1950) resonance 
was found. All subsequent fits were therefore tried with an S-wave threshold resonance 
plus MC shapes determined for contributions of a resonance at M = 2100 ± 20 (MeV/c 2 ), 
T = 160 ± 20 (MeV/c 2 ) (our optimum values), and phase space. The MC-determined event 
selection efficiency as a function of AM is shown in Fig. 20. The average efficiency, weighted 
by the threshold resonance contribution, as shown by the curves marked (1) in Fig. 21, was 
found to be e = 0.254. 
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FIG. 19. The AM = M(pp) — 2m p invariant mass distribution for the data for J/tjj 
curve illustrates the shape of the phase space contribution. 
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FIG. 20. MC-determined efficiency as a function of AM = M(pp) — 2m p for J/ip — > pp decays. 
The weighted average efficiency over the whole range is 25.4%. 



We find that if, like BES, we fit the AM distribution only in the region AM = 0-300 
MeV/c 2 , and do not consider the contribution due to the resonance at M(pp) = 2100 
MeV/c 2 , we obtain a good fit (x 2 /d.o.f. =15/26) which is essentially identical to that obtained 
by BES [5], with the results 

M(R thr ) = 1861^ 6 MeV/c 2 , 
r(i? thr ) = 0±f MeV/c 2 , 

B{J/4> -)• 7 i?thr) x B(R thT -)• pp) = (5.9+|;l) x (15) 

The errors are statistical only. 

Fitting data in the extended region AM = — 970 MeV/c 2 leads to a better determination 
of the phase space contribution; it is smaller than what is obtained if the fit is confined to the 
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FIG. 21. Fits of the AM = M{pp) — 2m p invariant mass distribution for J/tp — > *ypp decays. 
Dashed lines are the contributions of (1) i?thr> (2) /2(2100) and (3) the phase space. The solid line 
is the sum of all three contributions, (a) Fit in the full region AM = — 970 MeV/c 2 . (b) Same 
fit in the AM = -300 MeV/c 2 . 

small region, AM = — 300 MeV/c 2 . The fit with all parameters kept free (x 2 /d.o.f. =97/89) 
leads to the fractional contributions of /2(2100) of (22 ± 3)% and phase space (54 ± 3)%. 
The threshold resonance obtains 23llgg counts and leads to the following parameters 

M(R thr ) = 18371^1? MeV/c 2 , 
T(R thl ) = 0^ 4 MeV/c 2 , 

B(J/i; 7^) x B(R thv ->pp) = (11.41^1^) xlO- 5 . (16) 

The first errors are statistical, and the second errors are estimates of systematic errors 
obtained by varying the mass and width of the 2100 MeV/c 2 resonance by their la uncer- 
tainties. The same result, M(i? t hr) = 1837l^ 2 MeV/c 2 , is obtained when r(i2tkr) is fixed to 
20 MeV/c 2 . 

Figure 21(a) shows the fit in the region AM = — 970 MeV/c 2 , and Fig. 21(b) shows 
the same fit in the AM = - 300 MeV/c 2 region. 

The fit result for mass in Eq. 16 is consistent with the conjecture that the threshold 
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enhancement might be due to the tail of a subthreshold resonance at that mass. This 
possibility was raised earlier by BES with their observation of a resonance with mass M(R) = 
1833.7 ± 6.1 ± 2.7 MeV/c 2 in the reaction J/rj) -> *yR, R ->■ ir+ij-r] [18]. 



IX. SYSTEMATIC UNCERTAINTIES 

In Table VI we list various contributions to the systematic uncertainties in the branching 
fractions. References to previous CLEO studies for several of these are also given. The 
uncertainty in feed-down and continuum contributions to the background leads to system- 
atic uncertainties of 1%, 1.4%, and 2% in ^ypp, n pp, and rjpp, respectively. All the above 
contributions add in quadrature to 3.5%, 4.0%, and 5.5% uncertainty in •ypp, 7i°pp, and rjpp, 
respectively. Since the uncertainties in the fractions of individual resonance and continuum 
contributions are taken into account in the statistical uncertainties, the above are the sys- 
tematic uncertainties in the product branching fractions for the individual resonances in 
these decays. For B(ip(2S) — > 7pp), B(ip(2S) — > n°pp), and B(ip(2S) — > r/pp), the uncertain- 
ties in the fractions of individual resonance and continuum contributions lead to additional 
systematic uncertainties because the different contributions have different efficiencies. To 
take account of correlations in the fractions, the effective overall efficiencies were determined 
by MC simulations and the relative uncertainties found to be 2.4%, 0.8%, and 0.7% for ^pp, 
ir°pp, and r/pp, respectively. Thus the total systematic uncertainties in B(tp(2S) — > 7pp), 
B(i>(2S) ->■ n°pp), and B{^{2S) ->■ rjpp) are 4.2%, 4.1%, and 5.5%, respectively. 

The results for the branching fractions with systematic errors are given in Table VII. 

It is found that ±100 MeV/c 2 changes in the masses and widths of resonances introduce 
changes in branching fractions much less than 1%. 

TABLE VI. Sources of systematic uncertainties in branching fractions in %. References to previous 
CLEO publications for estimates of uncertainties are given in square brackets. 

Source [Ref] jpp ir°pp rjpp 

Number of ip{2S) [19] 2~!o 2X) 2.0 

Trigger efficiency [20, 21] 1.0 1.0 1.0 

Tracking efficiency [22] 2x1.0 2x1.0 2x1.0 



Particle identification [19] 


1.0 


1.0 


1.0 


7, 7T°, rj reconstruction [21] 


1.0 


2.0 


4.0 


Background subtraction 


1.0 


1.4 


2.0 


Sub total (quadrature) 


3.5 


4.0 


5.5 


Resonance fractions 


2.4 


0.8 


0.7 


Total (quadrature) 


4.2 


4.1 


5.5 



X. SUMMARY AND CONCLUSIONS 



Using CLEO data for 24.5 million ip(2S) we have studied decays ip{2S) — > 7pp, n°pp, 
and rjpp. In all three decays we find that intermediate N*(N*) states decaying into tt°N(N) 
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TABLE VII. Summary of the measured quantities. First errors are statistical, second errors are 
systematic. 



Quantity 






Events 


Our result 
(10- 5 ) 


PDG08 

(10- 5 ) 


o{ip{zb) - 


> jpp) 




348 ±22 4.18 ±0.26 ±0.18 


z.y it u.o 


oyipyzs) - 






948 ± 37 


15.4 ±0.6 ±0.6 


1004.1 7 

lo.o It 1.1 


&{1p{ZS) - 


* wpp) 




154 ± 16 


5.6 ±0.6 ±0.3 


(5(14-1 O 
D.U It l.Z 


Bty(2S) - 


> 7/2(1950)) x B(/ 2 (1950) -> 


pp) 


111 ± 19 


1.2 ±0.2 ±0.1 




t5{ip{Zb) - 


i. f foi v Rf f (O^ i. 

* 7J2i^it>UJJ X D(j2(ZLbU) — > 


pp) 


73 ± 18 


0.72 ±0.18 ±0.03 




Bms) - 


> piVf(1440)) x B(JVf(1440) 


->■ P7T°) 


474 ± 42 


8.1 ±0.7 ±0.3 




Bty(2S) - 


* piV 2 *(2300)) x £(iV 2 *(2300) 


^pvr ) 


265 ± 39 


4.0 ±0.6 ±0.2 




B{^{2S) - 


> 7T #l(2100)) X £(#1(2100) 


— >■ pp) 


76 ±29 


1.1 ±0.4 ±0.1 




Bty(2S) - 


» ^°i? 2 (2900)) x £(#2(2900) 


^pp) 


133 ± 38 


2.3 ±0.7 ±0.1 




BtypS) - 


* piV*(1535)) x £(iV*(1535) 


-¥ prj) 


123 ± 16 


4.4 ±0.6 ±0.3 




Bms) - 


* ^1(2100)) X £(i?i(2100) - 


^pp) 


31 ±10 


1.2 ±0.4 ±0.1 





and rjN(N), and fj, aj meson resonances decaying into pp make important contributions 
to the total decay. We have determined branching fractions for the total decay and for 
the contributions of the individual intermediate states. For the total decays our branching 
fractions have factors two to three smaller uncertainties than in the current literature. The 
product branching fractions for decays through individual intermediate states have been 
determined for the first time. The results are summarized in Table VII. 

We do not find any evidence for a threshold enhancement in any of the three ip{2S) decay 
channels. For ip(2S) — > •ypp we set a stringent upper limit for the threshold resonance -R t hr, 
B{^{2S) ->■ 7#thr) x B{R thT ->■ pp) < 1.6 x 10- 6 at 90% CL. 

With a limited sample of 8.6 million J/ip available to us from ip(2S) — > n + iT~J/'ip 
we have searched for J/ip — > 7-R t hr- We find a pp threshold enhancement. When it 
is analyzed taking into account an enhancement at M(pp) = 2100 MeV/c 2 , we obtain 
M(R thr ) = 1837+1°1| MeV/c 2 , T(R thr ) = 0^ 4 MeV/c 2 , and B{J/^ ->■ 1 R thr ) x B(R thr -> 
PP) = xlO- 5 . 
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